Discharge flow of granular particles through an orifice on a horizontal hopper: Effect of the hopper angle
Wang Xin, Zhu Hong-Wei, Shi Qing-Fan, Zheng Ning
School of Physics, Beijing Institute of Technology, Beijing 100081, China

 

† Corresponding author. E-mail: qfshi123@bit.edu.cn ningzheng@bit.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11475018 and 11974044) and the National Key Research and Development Program of China (Grant 2016YFC1401001).

Abstract

We experimentally investigate the effect of the hopper angle on the flow rate of grains discharged from a two-dimensional horizontal hopper on a conveyor belt. The flow rate grows with the hopper angle, and finally reaches a plateau. The curve feature appears to be similar for different orifice widths and conveyor belt-driven velocities. On the basis of an empirical law of flow rate for a flat-bottom hopper, we propose a modified equation to describe the relation between the flow rate and hopper angle, which is in a good agreement with the experimental results.

PACS: ;45.70.Mg;;47.56.+r;
1. Introduction

The improvement in flow rate of discrete particles passing through a bottleneck is of importance for plenty of fields, widely ranging from industrial transportation, traffic flow to evacuation of panicking crowds.[13] Usually, controlling the orifice size is a direct and convenient way to obtain the desired flow rate. However, when the orifice size cannot be changed in some particular situations, the optimization of flow rate under this constraint is a practical issue worthy of study.

In a vertical gravity-driven hopper flow, some methods have been adopted to improve the flow rate. Obviously, widening the orifice is the most effective way. It is found that the packing fraction and velocity of particles near the orifice simultaneously increase with the orifice size.[4] Changing the exit position can also increase the flow rate. Numerical results show that the flow rate remains constant when the exit is far from the wall and increases exponentially when the exit moves close to the wall.[5] Shaking the hopper is another method frequently used to enhance the flow rate.[6,7] By an external vibration, the arches near the orifice that block the flow become unstable and are readily broken. Also, an obstacle placed above an orifice with an optimal distance is used to suppress the formation of arches.[812] When the obstacle size and position are properly chosen, the flow rate can be increased up to 10%. The hopper angle has been proved to be another effective control parameter for the enhancement of the flow rate.[1317] Based on an hourglass theory, a formula was derived for the flow rate of a narrow-angled hopper.[1] The formula gives an excellent qualitative prediction of the measured dependence of the mass flow rate on the density, orifice diameter, hopper angle, and angle of internal friction. However, it does not provide an accurate quantitative prediction, and the theoretical value from the formula is about twice the experimental measurement. Two of the reasons for the discrepancy are the assumptions that the hopper walls are smooth, and that the stresses fall to zero on a free-fall arch. The concept of free-fall arch is not proved by experiment and this could lead to unquantifiable uncertainties in the prediction. Later, the flow rate from an angled hopper with an obstacle above the orifice was further investigated on the basis of the researches above. The optimal flow rate was increased by up to 16% compared with that without an obstacle. Using the hourglass theory and a velocity-density relation, an empirical law that relates the flow rate, hopper angle and the obstacle position was proposed.[19]

For a horizontal hopper flow, the previous work mainly focused on a flat-bottom hopper. The orifice size, the driven velocity, and the initial packing fraction of particles have been considered as key factors for the improvement of the flow rate.[2022] The effect of obstacles has been explored, showing that the obstacles on both sides of an orifice can effectively increase the flow rate by up to 15.6%.[24] In contrast to the intensive studies on the effect of hopper angles in the vertical gravity-driven system, the flow rate dependence on the hopper angle in a horizontal system still remains unknown, although this type of hopper flow, such as that driven by a conveyor belt, has wide applications in industry. In addition, traffic jam[25] or emergency evacuation management[3] could even draw inspiration from the optimal design of the orifice geometry.

In this work, our aim is to present a parametric study of a granular flow driven by a conveyor belt through a horizontal silo with different hopper angles. We experimentally report the effect of the hopper angle on the flow rate. The flow rate dependence is measured with different orifice sizes and conveyor belt-driven velocities. On the basis of an empirical law of flow rate from a flat-bottom vertical hopper, a modified law is proposed to directly relate the flow rate to the hopper angles. The potential application for the unjamming of particles near an orifice by means of the hopper angle is briefly discussed.

2. Experimental details

The experimental establishment consists of a conveyor belt above which a two-dimensional, hourglass-shaped frame is fixed, with 1750 Plexiglas disk particles confined in this frame, as shown in Fig. 1. A motor drives the belt at a constant velocity ranging from 3 cm/s to 15 cm/s. A rectangular orifice with width D is located at the center of the bottom wall of the frame. The hopper angle between the horizontal and the hopper bottom can be varied from θ = 0° (flat-bottom hopper) to θ = 90° (pipeline hopper). The initial packing fraction of the dense-packing disks is set to be ϕ0 = 0.86 ± 0.02 by using a previous protocol.[24] A steady flow status is always maintained since the orifice is larger than a critical size at which clogging occurs (D ⩾ 6d, d is the diameter of the disk particle).[23]

Fig. 1. (a) Schematic diagram of the experimental establishment. Disk particles on a conveyor belt are driven to flow through an orifice. A high-resolution camera is used to record the number of the particles in the frame. (b) Illustration of the hopper angle. By adjusting the position of the bottom wall (highlighted by white lines), different hopper angles from 0° to 80° can be obtained. d = 11 mm is the diameter of the disk particle used in the experiment.

A high-resolution video camera hung above the frame records the temporal evolution of the disk number inside the frame, and thus an average flow rate during the steady flow is extracted as before.[20,24] The measurement is repeated 6 times for each data point, and the flow rate is computed as an arithmetic average of these trials.

3. Results and discussion

We report that the dependence of the flow rate on the hopper angle with a given orifice width D = 10d, as shown in Fig. 2(a). Here we examine the dependence curve for a wide range of belt-driven velocities; these flow rate–hopper angle curves look similar and the curve feature is robust: the flow rate monotonically increases with the hopper angle, until it reaches a plateau. In order to show the dependent feature more clearly, the data in Fig. 2(a) are normalized by the belt-driven velocity (Fig. 2(b)). Compared with the flow rate from a flat-bottom hopper (θ = 0°), the flow rate increases by about 13%–15% when θ > 50°. It is worth pointing out that the improvement in the flow rate solely due to the geometrical modification of the orifice is not a marginal effect. This is because the theoretical maximal value of the flow rate is ,[20,23] where ϕ is the initial packing fraction of the dense-packing particles (note ϕ is not the packing fraction of flowing particles near the orifice), V the belt velocity, D the orifice width, and d the particle diameter. If the orifice size and the belt velocity are fixed, the theoretically ideal flow rate is just 23% larger than that for a flat-bottom hopper, where the flow rate of a flat-bottom hopper ,[20,23] k = 1.86 for our measurement. Generally, the value of k was in a range of 1–3.5 in previous literature. The value of k has been found to be dependent on the particle and hopper properties in a range of 1 < k < 2.[26] If the flowing particles consist of sands, the value of k can reach 2.9.[27] For a horizontal flat-bottom hopper flow system, k = 2.1 ± 0.2 is similar to that in our system.[20] If the conveyer belt velocity is very large, the flow rate is related to the friction between the conveyer belt and flowing particles, and k = 3.23.[28]

Fig. 2. (a) The flow rate Q versus hopper angle θ for different conveyor belt velocities. The orifice width D = 10d. The size of the error bars is almost the same as that of the symbols. (b) Dependence of normalized flow rate on the hopper angle shows the curve feature in (a) more clearly.

In the previous literatures, an empirical law equivalent to the 2D Beverloo’s law,[23] , was proposed to describe the flow rate of particles on a conveyor belt passing through an orifice.[20,23,24] However, the law is only applicable to a flat-bottom horizontal hopper on a conveyor belt. To interpret the results in Fig. 2, we adopt a modified model for a two-dimensional flat-bottom vertical hopper proposed by Mancok and Janda et al.[4,29] The flow rate Q in this model is derived based on the self-similar velocity and density profiles at the orifice, , where C′ is a constant depending on the particle diameter and on the curvature of the density profile, R is the half-width of the orifice, ϕ0 is the asymptotic value of the packing fraction for big orifices, α and β are fitting parameters. Due to the similar self-similarity in the velocity and density profiles at the orifice for the horizontal hopper (see the supplementary materials), we assume that the flow rate in the two-dimensional flat-bottom horizontal hopper shares the same mathematical form as that in the vertical hopper,

where C′, ϕ0, α, and β have the same physical meanings as those in the vertical hopper model, D = 2R is the orifice width, and vb is the belt velocity. To examine the speculation, the reduced flow rate is plotted as a function of the orifice width in units of the particle diameter d, (d = 1.1 cm is a constant in our experiment), as shown in Fig. 3. The solid line in this figure is the fitting curve with an exponential form: C′(1–αeβD/d), where C′ = 0.95, α = 0.74, and β = 0.19. An excellent agreement between the experimental flow rate and the prediction in Eq. (1) confirms the assumption that the flow rate can be well described by Eq. (1).

Fig. 3. Reduced flow rate Q* as a function of the orifice width D in units of the particle diameter d. The solid line is a fitting curve which follows the form of an exponential function C′(1–α eβD/d), where C′ = 0.95, α = 0.74, and β = 0.19.

On the basis of the modified flow rate expression for the two-dimensional flat-bottom vertical hopper , a physical model to describe the effect of the hopper angle on the flow rate has been proposed.[16] In this model, the hopper angle is considered to be influential in the radius of curvature of an arch at the orifice. The new flow rate expression is , where d is the particle diameter and θ is the hopper angle. Due to the similar geometrical feature at the orifice between the horizontal and vertical systems, we assume that the mathematical form of this law still holds for the horizontal hopper,

In Fig. 4, the experimental flow rates for different hopper angles and orifice sizes are scaled to compare with the prediction of Eq. (2). The data points collapse fairly well and reveal an excellent agreement with the fitting curve from Eq. (2), which suggests the expression in Eq. (2) can describe the flow rate well for all hopper angles.

Fig. 4. Scaled flow rate as a function of D/(d cos θ) for different belt velocities and orifice sizes. The solid symbols are for an orifice width D = 6d = 6.6 cm, and the hollow symbols for D = 10d = 11 cm. All data points collapse together to form a master curve. The solid curve is a fitting line which follows an exponential function Cx(1–α eβx), where , C′ = 0.95, α = 0.74, and β = 0.19.

Recent work pointed out that the hopper angle may significantly affect the clogging probability of a particle through an orifice. By increasing the hopper angle, the clogging probability can be reduced by three orders of magnitude.[26] Therefore, the increase of the hopper angle leads to a resultant decrease of temporary and permanent clogging arches which consist of a few clogging particles developing at the orifice. An appropriate increase of the hopper angle will probably maintain the flowing status of particles without a complete arrest, even if the orifice size is smaller than the clogging size in a flat-bottom hopper.

4. Conclusion

We present experimental results of the effect of the hopper angle on the flow rate of grains discharged from a two-dimensional horizontal silo on a conveyor belt. The flow rate grows with the hopper angle, and finally reaches a plateau. Based on an original empirical law, we propose a modified equation which can well describe the dependence of the flow rate on the hopper angle within a wide parameter range. The results reported in this work show a practical method for the improvement of the flow rate through a bottleneck, especially for the case where the orifice size and driven velocity have to be restricted. In addition, an increase in the hopper angle will probably be beneficial to maintain flowing of discrete objects through a narrow exit even when the bottleneck is smaller than a critical size of clogging in a flat-bottom hopper. The feature is probably not only applicable to the conveyor belt system, but to other systems of discrete particles flowing through a constriction such as pedestrian, active agents, and mirco- or nano-particles.

Reference
[1] Nedderman R M 1992 Statics and Kinematics of Granular Materials Cambridge Cambridge University Press 301 10.1017/CBO9780511600043
[2] Helbing D Johansson A Mathiesen J Jensen M H Hansen A 2006 Phys. Rev. Lett. 97 168001
[3] Helbing I F D Vicsek T 2000 Nature 407 487
[4] Janda A Zuriguel I Maza D 2012 Phys. Rev. Lett. 108 248001
[5] Wang T W Li X Wu Q Q Jiao T F Liu X Y Sun M Hu F L Huang D C 2018 Chin. Phys. 27 124704
[6] Chen K Stone M B Barry R Lohr M McConville W Klein K Sheu B L Morss A J Scheidemantel T Schiffer P 2006 Phys. Rev. 74 011306
[7] Wen P Zheng N Nian J Li L S Shi Q F 2015 Sci. Rep. 5 09880
[8] Zuriguel I Janda A Garcimartín A Lozano C Arévalo R Maza D 2011 Phys. Rev. Lett. 107 278001
[9] Lozano C Janda A Garcimartín A Maza D Zuriguel I 2012 Phys. Rev. 86 031306
[10] Yanagisawa D Kimura A Tomoeda A Nishi R Suma Y Ohtsuka K Nishinari K 2009 Phys. Rev. 80 036110
[11] Zuriguel I Olivares J Pastor J M MartínG ómez C Ferrer L M Ramos J J Garcimartín A 2016 Phys. Rev. 94 032302
[12] Yang S C Hsiau S S 2001 Powder Technol. 120 244
[13] Anand A Curtis J S Wassgren C R Hancock B C Ketterhagen W R 2008 Chem. Eng. Sci. 63 5821
[14] Choi J Kudrolli A Bazant M Z 2005 J. Phys: Condens. Matter. 17 S2533
[15] Uñac R O Vidales A M Benegas O A Ippolito I 2012 Powder Technol. 225 013
[16] Benyamine M Aussillous P Dalloz D B 2017 EPJ. Web Conferences 140 03043
[17] Yang L N Chen Y Q 2019 Chin. Phys. Lett. 36 024501
[18] Wan J Wang F Yang G Zhang S Wang M Lin P Yang L 2018 Powder Technol. 335 147
[19] Alonso-Marroquin F Azeezullah S I Galindo-Torres S A Olsen-Kettle L M 2012 Phys. Rev. 85 020301
[20] Aguirre M A Grande J G Calvo A Pugnaloni L A Géminard J C 2010 Phys. Rev. Lett. 104 238002
[21] De S B Xun S Z Guang L X Zheng Q P Xiao W T Kun Q L 2003 Phys. Rev. 67 062301
[22] Cordero M J Pugnaloni L A 2015 Powder Technol. 272 290
[23] Aguirre M A De Schant R Géminard J C 2014 Phys. Rev. 90 012203
[24] Zhu H W Wang L P Shi Q F Li L S Zheng N 2019 Powder Technol. 345 053
[25] Helbing D 2001 Rev. Mod. Phys. 73 1067
[26] Nedderman R M Laohakul C 1980 Powder Technol. 25 91
[27] Fowler R T Glastonbury J R 1959 Chem. Eng. Sci. 10 150
[28] Manuel J C Luis A P 2015 Powder. Technol. 272 290
[29] Mankoc C Janda A Arévalo R Pastor J M Zuriguel I Garcimartín A Maza D 2007 Granul. Matter. 9 407